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The influence of electronic structure evolution upon pressure on the temperature dependencies of 
electrical resistivity of pure Np, Pu, Am, and Cm metals have been investigated within coherent 
potential approximation (CPA) for many-bands conductivity model. Electronic structure of pure 
actinide metals was calculated within the local density approximation with the Hubbard U and 
spin-orbit coupling corrections (LDA-f (/+SO) method in various phases at normal conditions and 
under pressure. They were compared with the corresponding cubic (bcc or fee) phases. The densities 
of states of the latter were used as a starting point of model investigations of electrical resistivity. 
The obtained results were found in good agreement with available experimental data. The nature 
of large magnitude of resistivity of actinides was discussed in terms of the proposed conductivity 
model using the ab initio calculated parameters. 

PACS numbers: 71.27.-(-a, 71.10.-d, 71.20.-b 



I. INTRODUCTION 

For many decades transuranium metals have been of 
interest to scientific community due to the unique combi- 
nation of structural, electronic, and kinetic propertiesii 
Their rich phase diagrams^ include a number of struc- 
tural transitions upon pressure or increasing the temper- 
aturei^^— Intermediate or jj coupling schemes, in con- 
trast to usual Russel-Saunders [LS) one, were found to 
be more appropriate for 5/ electrons in transuranium 
metals due to the presence of large spin-orbit coupling 
in actinidesi^"— Additionally, magnetic stabilization of 
phases was found in curiumi^iii Drastic changes of crys- 
tal and electronic structures and their interplajJ^ do not 
deplete this list of unusual properties of the metals. 

Anomalous resistance properties in transuranium met- 
als are also well knownJ^ The negative temperature co- 
efficient of resistivity (TCR) was measured in a- and 
i5-Pu and some Pu-based alloys at high temperatures 
(T > 6'£))i^ii^ Moreover, p ^ T^ dependence was found 
in these compounds at low temperature F^^ii^iii Electri- 
cal resistivity (ER) behavior in pure Np and Am met- 
als is ordinaryfi^^— namely, p ~ T in high and p ^ T^ 
and p ~ 2^4-5 jj^ j^.^ temperature ranges, respectively)^ 
The superconducting state was found in Am at ambient 
pressure with the characteristic temperature Tc ~ 0.5 - 
2.0 K32, Curium metal also demonstrates ordinary be- 
havior of ER typical for antiferromagnetic metal with 
p ~ T^ at low temperature and p ~ T at high temper- 
ature, respectively. However, temperature dependencies 
of ER in transuranium metals show anomalously high 
values of ER (typically, 80 -;- 140 fikVi ■ cm) at room tem- 
perature. For comparison, ER of transition metals at the 
same conditions are about 1 -;- 10 p,k^ ■ cm. 

Such high resistivity values of trasuranium metals have 



no reasonable explanation yet and all previous model in- 
vestigations substantially underestimate the experimen- 
tal values of ER. For example, Mott two-band conduc- 
tivity model combined with coherent potential approx- 
imation (CPA) shows, that Np, Pu, Am, and Cm are 
typical metals at high temperatures,^'^ and observed ER 
behavior is described well in relative units. However, ab- 
solute values of ER in transuranium metals cannot be 
estimated from these calculations. Indeed, the estima- 
tion of maximum value of ER - taking into account only 
interband transitions of scattering conductivity electrons 
with additional accounting that probabilities are propor- 
tional to the DOSes values at the Fermi level - gives the 
ER values at most 20 -;- 40 pkQ. ■ cm, that underestimates 
the experimental data by a factor of 3. 

It is well known, that the high-resistivity values can 
be provided by specific Kondo-like resonance. Actinides, 
usually considered as heavy-fermion systemSfi^ should 
demonstrate appropriate kinetic and magnetic proper- 
ties in low temperature range. However, we have no 
experimental data for the detailed discussion in terms 
of such a model. Previously reported investigations of 
ER within common Kondo mode l^^'^^ have a number of 
problems and ambiguities^^ due to selecting of the "mag- 
netic" part of ER. A principal problem arose if one uses 
spin-fluctuation modeP- for describing ER of (5-Pu and 
other actinides. Indeed, one can find good agreement 
of calculated and experimental data only adjecting the 
Stoner factor..^^ However, such a high magnitude of the 
Stoner factor should result in specific temperature behav- 
ior of magnetic susceptibility, typical for amplified para- 
magnets. In contrast to these expectations, ordinary for 
paramagnets temperature dependence of magnetic sus- 
ceptibility was measured in the systems under consider- 
ation!^ 

On the other hand, recent NMR experiments for 5- 



Pu revealed a combination of Curie- Weiss and van-Fleck 
temperature behavior of spin susceptibility. ^!^-' That 
demonstrates strong influence of spin-density fluctua- 
tions on electronic structure and can provide anoma- 
lies of electron heat capacity in the same temperature 
range, as reported in Ref. |30- Actually, strong spin- 
density fluctuations effects can play a noticeable role in 
actinides, providing p ~ T^ in very low temperature re- 
gion, and suppressing superconductivity. Note also, that 
the p '^ T"^ dependence and negative TCR (observed in 5- 
Pu) can be described in terms of strong electron-phonon 
coupling j^^i^^ Thus, theoretical explanation of the high 
values of resistivity of transuranium metals will require 
significant efforts, and hence ab-initio calculations are of 
high interest for the understanding of ground state prop- 
erties in these materials. 

Recent experiments and electronic structure calcula- 
tions revealed numerous effects due to strong correla- 
tions of the 5/ electrons in Np, Pu, Am, and CmM 
A number of band methods and approximations have 
been applied to describe magnetic and spectral proper- 
ties of transuranium metals ji^ Nonmagnetic ground state 
of pure plutonium metal observed experimentally^^, was 
reproduced in the electronic structure calculations by the 
LDA-I-C/-I-SO methodi^ There the exchange interaction 
was found to be the reason of artificial antiferromagnetic 
ordering in previous LSDA-I-C/ investigations. Also non- 
magnetic ground state was obtained in around-mean-fleld 
version of the LBA+U method^^ then in LDA -t- Hub- 
bard I approximatiou )'^^''^^ and hybrid density function- 
als with a dominant contribution of HF functionali^ Re- 
cently, the reliability of these results was proved by the 
more detailed analysis of exchange interactiou j'^^i'^^ and 
also by the LDA+DMFT method,'"' combining the LDA 
approximation with the Dynamical Mean-Field Theory 
(DMFT) in various modifications.—"— However, consis- 
tent interpretation of spectroscopic data by these calcu- 



lations is not found yet 
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In this paper we report the results of electronic struc- 
ture and resistivity calculations for transuranium met- 
als at normal conditions and under pressure. In terms 
of CPA the many-band conductivity model was derived 
and applied for fcc-Pu, Am, and Cm and bcc-Np numer- 
ical investigations of electrical resistivity using DOSes of 
metals calculated within the LDA-|-[/-|-SO method as a 
starting point. The derived model allows one to account 
for initial DOSes modifications in a direct way as a re- 
sult of temperature erosion and its renormalization due 
to interband s — >■ d, s — >■ /, d — >■ /, and / — >■ d electron 
transitions. The DOSes evolution of the metals with tem- 
perature at normal conditions and under pressure were 
presented and discussed. Results of the ER calculations 
are found in good agreement with available experimen- 
tal data. In terms of the proposed conductivity model, 
using ab initio obtained parameters, the nature of high- 
resistivity values in actinides is also discussed. 
II. ELECTRONIC STRUCTURE 
CALCULATIONS 



We investigated the electronic structure of transura- 
nium metals within the LDA+C/+SO method described 
in detail in Ref. 33;. In this method the exchange interac- 
tion (spin polarization) term in the Hamiltonian is imple- 
mented in the general nondiagonal matrix form regard- 
ing the spin variables. This form is necessary for correct 
description of 5/ electrons for the case of jj and inter- 
mediate couplings. ^'^ The LDA-f /7-l-SO method is intro- 
duced by the effective single-particle Hamiltonian which 
supplements initial local density approximation (LDA) 
functional by a term with orbital dependent potential 
V^^, , which accounts for Coulomb interactions (^ and 



spin-orbit coupling termi^ 
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Here U and Jh are screened Coulomb and Hund ex- 
change parameters which are determined in the constrain 
LDA calculations j^2i^ The screened Coulomb interac- 
tion matrix elements {m,m" \ Vee \ m' ,m"') could be 
expressed via these parameters4^ In the constrain LDA 
procedure a screened Coulomb interaction of 5/ electrons 
is evaluated that requires the choice of screening chan- 



nels taken into account in the constrain LDA calcula- 
tions. Taking s, p, and d channels, for Np, Pu, Am, and 
Cm we obtained the Coulomb parameter value [/ « 4 eV 
in good agreement with the previously used value.— 

In contrast to the direct Coulomb parameter f/, the 
exchange parameter Ju is evaluated as a difference of in- 
teraction energy for the electrons pairs with the opposite 



and the same spin directions. Here the screening process 
is defined by the charge but not spin state of the ion, 
then the screening contribution is canceled for exchange 
Coulomb interaction and parameter Ju does not depend 
on the choice of screening channels. For neptunium and 
plutonium, the value of Hund exchange parameter Jh 
was reported to be Jh = 0.48 eV,— for americium - 
0.49 eV, and for curium - 0.52 eV,^i^ 

In Eq. (m the off-diagonal exchange interaction terms 
with n^f,^/,/ are significant for actinide elements. It was 
demonstrated in Ref. 33, that omission of these terms can 
result in incorrect antiferromagnetic ground state for fcc- 
plutonium metal. This fact is due to the intermediate or 
jj coupling scheme taking place in most actinide metals, 
including Pu, Np, and Cm. But in pure Am metal jj 
coupling is present. It means that total moment J is 
well defined, but not spin S and orbital L moments, as 
in usual LS coupling scheme. In this case, the basis of 
eigenfunctions of total moment operator {jruj} is the 
best choice. The matrix of spin-orbit coupling operator 
is diagonal in this basis but not the exchange interaction 
(spin-polarization) term in the Hamiltonian. 

Whereas in the LS coupling scheme S and L operators 
are well defined. Then the basis of LS orbitals, which 
are eigenfunctions of both spin S and orbital moment L 
operators, is a good choice. In this case it is possible to 
define quantization axis in the direction of spin moment 
vector so that occupation and potential matrices will be 
diagonal in spin variables. 

When the intermediate coupling is realized, neither jj 
basis, nor LS is valid, and occupation matrix is nondiag- 
onal in both orbital bases and both terms in the Hamilto- 
nian: spin-orbit coupling and exchange interaction, must 
be taken in a general nondiagonal matrix form, and some 
finer treatment is necessary. Intermediate basis is still 
cumbersome and is used for model atomic calculations. 

For the actinide metals under consideration the 
LDA-(-C/-|-SO densities of states (DOS) were taken as a 
starting point for further CPA simulations of ER temper- 
ature dependencies within the model described below. 



III. COHERENT POTENTIAL 

APPROXIMATION FOR MANY BAND 

CONDUCTIVITY MODEL 

In the general case, the many-hand conductivity model 
is suitable for explanation of kinetic properties in ac- 
tinides due to equal probabilities of transitions of (s -\-p) 
electrons into almost empty d and / bands, since the val- 
ues of DOSes at the Fermi level in d and / bands are 
close to each other. Thus, ratios between partial d and 
/ DOSes at the Fermi level lead to opening of additional 
channel of direct d ^ f and back f ^ d electrons transi- 
tion. Moreover, the evolution of the d and / DOSes with 
temperature and direct and back transitions depends on 
the evolution of DOSes of the other shell. Then it is 
reasonable to obtain DOSes at finite temperatures with 



regular self-consistent procedure and derive correspond- 
ing CPA set of equations without model simplification, in 
the same way as it was made previously in Ref. |3l|, when 
only direct s — >■ d and s — > / transitions were accounted 
for. 

We start considering the s{p), d, and / electrons per- 
forming intra- and inter-band transitions as a result of 
their scattering at long wave phonons. We assume also 
that the accepting d or / bands are partially filled. Then 
the Hamiltonian of electron subsystem H = Ho + Hint 
can be written down in the following form: 
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H = ^i?;a+a, + ^ E e-'^'-''-''-^Vu'{n)d+di,, (3) 

where Ei is a periodical part of electrons energy. Com- 
bined index I includes the band index j [j ~ s, d, /) and 
wave vector fc; _R„ is a radius- vector of the n-th site of a 
crystal lattice. Operator Vij'{n) describes the intensity 
of electron-phonon interaction. If thermal displacements 
of ions are small, the operator Vi^i'{n) can be written as 



where 
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is the parameter of intensity of intra- and inter-band 
transitions due to electron-phonon scattering. M and 
S are the mass and the sound velocity in metal, respec- 
tively; go is the maximum value of q and Kp is the Fermi 
wave number of electron; A;.// are the (Bloch) parameters 
of electron-phonon coupling. 

Using ab-initio obtained DOSes of metal as a starting 
point of numerical calculations we assume that effects of 
s, p, d, f hybridization are accounted for in the elec- 
tron ground state. For simplicity we keep the same band 
notations after renormalization. 

A set of CPA equations can easily be derived using 
Dyson equation and definition of T-matrix. Let us de- 
termine the total resolvent of the full energy operator 
H: 



i?= (z-H)- 



(6) 



and the strictly diagonal part of the total resolvent R in 
the Hq representation: 



G=iz-Ho- A)- 



(7) 



where A is a strictly diagonal shift operator in the Hq 
representation. Broadening of single electron levels is 
described as: 



1 



A^Ay 



i{k-k',R„) 



^A 
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(8) 
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The real part of the coherent potential Aj determines 
the shift rjj, and its imaginary part 7 determines the 
broadening of single-electron levels. 

To derive equations within many-band CPA, the Dyson 
identity was used: 



R^G + G{V~'A)R. 



(9) 



Scattering operator T can be determined in a convenient 
form as: 



brackets should be taken. The series p2)) contains mu- 
tually compensated block terms, in which the shift oper- 
ator is included in indirect form.— Excluding the com- 
pensated block terms from Eq. (J12p . and averaging over 
ion displacements, one finally obtains the shift operator 
as 



A 



VGV + VGVGV 



(13) 



R — O — CtTCj. 



(10) 



Multiplying both parts of identity p^ by G ^ , one ob- 
tains the following expression for the scattering operator: 



f = G-^{R-G)G-^. 



(11) 



Using the Dyson identity © and Eq. pT|) for the shift 
operator, the following operator series can be written: 



A = [{V - l\)G{V - k) + 
{V - k)G{V - A)G(\/ - A) + . . .Uag- 



(12) 



Here the diag index means that the diagonal part in the 
Hq representation of the sum of operator products in 



Here only strictly diagonal terms in the Hq representa- 
tion ([...]£)) are accounted for, and items containing the 
blocks are omitted. The brackets < ... > mean aver- 
aging over phonons. Let us assume also, that at high 
temperature the operator Vw {u) can be replaced with 
the one averaged over wave vectors k and k' functions of 
fluctuating variable u ~ Vn.jj'{u)^ Then the series P^ 
can be summed up accurately (in the convergence range 
|Vj(w)Fj| < 1) as it is shown in the Appendix for single- 
electron and single-site approaches. Using for simplifica- 
tion also \Vn,sj{u)Fs\ < 1, \Vn,sFs\ / 1 14,(d)/^(d)/ | < 1, 

/\VjjiFj\ <^ 1 and notation Vnjj{u) = U„j(u), 



WsjF., 



for the averaging over ion thermal displacements s- 
d-band coherent potentials, one obtains: 



and 



J 



E 



duP{u 



V^{u)F^[l - Vr{u)F^,] + Vsd{u)Vsf{u)Vdf{u)F.^F^, 
[1 - Vd{u)F,][l - Vf{u)F^] - V;if{u)F,F^ 



(14) 



A„ 



duP{u 



Vi{u)F^[l ~ Vf(,u)Ff] + V^f{u)Ff[l + Vd{u)F^] 



[^-Vr{u)F,][l-Vf{u)Ff 



ylM)F,Fj 



(15) 
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-u'^/2P 



(3 = Zo.^jj,2T/9d (16) 



is the Gauss distribution functionS^i^ and 9d is the De- 
bye temperature. Equation for coherent potential of the 
/ band electrons is the same as Eq. ([T5|) but a replace- 
ment of band indices s <=^ / is necessary. 

Note also, that an assumption \Vn,sj{u)Fs\ ^ 1 in 
Eqs. (|14l) and (|15l) means physically that the partial fill- 
ing of (i(/)-bands sets the traps for mobile s-electrons 
and corresponds to the Mott idea, successfully applied 
previously for two band conductivity model. 

On the other hand, Eq. ([T5|) describes direct and back 
transitions of d(/)-electrons and provides a possibility of 
accounting for d(/)-bands modifications with tempera- 
ture. These equations are solved within self-consistent 
loop. A solution of these equations allows one to calcu- 
late coherent potential for s-conductivity electrons. 



Simple assumption leads to the previously obtained re- 
sults. Indeed, using for example Vjji — 0, i.e., neglecting 
the inter-band transitions, three independent equations 
corresponding to the single-band model of the CPA can 
be obtained from Eqs. (fTl)) and (fT5|) ^ 



+00 



(A,) = 



du P{u) 



V,Hu)F, 
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Also, one can found the result for two band conductivity 
model and perturbation theory series for the coherent 
potential. 

Accounting for this renormalization, and also for the 
values of DOSes of bands and their modification with the 
temperature, one can estimate the magnitude of ER in 
actinides. Note, that effects of irradiation will not be 
taken into account in this result. 
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FIG. 1: (Color online) Partial 5/ densities of states for a- and /3-Np (left panels) and a- and 7-Pu (right panels) in real crystal 
structures are compared with the bcc-(Np) and fcc-(Pu) DOSes for the corresponding cell volumes per actinide atom (see 
Ref . ISa) . The Fermi energy corresponds to zero. 

TABLE I: Symmetry groups of Np, Pu, Am, and Cm and results of the electronic structure calculations for the real crystal 
structures within the LDA+f/+SO method. The largest values of occupation matrices off-diagonal elements (OD) in {1/5} and 
{j'mj} basis sets are given (see details in Ref. 1331). Then the calculated values for spin (5), orbital (L), total (J) moments, 
Lande factor, and effective magnetic moment (/il'f'f) in /is are presented. ^^ 
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6.42 



IV. RESULTS 

Ground state of metals at normal conditions 
and under pressure 



LDA+t/+SO method. As one can see from the anal- 
ysis of ofF-diagonal occupation matrix elements, in all 



Electronic structure of Np, Pu, Am, and Cm in real 
crystal structures (Table 1) were calculated within the 
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FIG. 2: (Color online) Partial 5/ densities of states for Ami, AmIII, and Am IV (left panels) and Cml, Cmlll, and Cm IV 
(right panels) in real crystal structures are compared with the fcc-DOS for the corresponding cell volumes per actinide atom 
(see Ref . ISa) . The Fermi energy corresponds to zero. 



metals (except S-Pu and Ainll) in both LS and jrrij 
bases off-diagonal elements are substantial and compara- 
ble. That means that intermediate coupling takes place 
in these metals. Simple cubic S-Pn and Amll are well de- 
scribed within the jj coupling scheme that is confirmed 
by /^ electronic configuration and smaller hybridization 
of j = 5/2 subband with j = 7/2 subband and reflected 
in its smaller occupation. 

In all our calculations for Pu and Am only nonmagnetic 
solution was found, whereas in neptunium effective mag- 
netic moment was obtained about 3 /is and depends on 
phase. In curium effective magnetic moment strongly dif- 
fers in the phases, see Table 1. Experimental value of the 
moment in Cm was measured as 7.85 ^b^^^ While model 
calculations predict magnetic moment from 7.94 fis in 
ionic picture in the assumption of pure LS coupling to 
7.6 fiB in the assumption of intermediate coupling,Sl 

For the many-band conductivity model presented in 
Section IIIII partial densities of states from ab initio cal- 
culations are used. In Figs. [1] and [2] we present 5/ partial 
densities of states for the calculated crystal phases. Re- 
sults for cubic (7-Np, S-Pn, Amll, and Cmll) structures 
under pressure were reported elsewhere.~ In all DOSes 
one can distinguish two groups of peaks attributed to 
the subbands with the total moment value j — 5/2 at 
the lower energies and j = 7/2 at higher ones splitted 



by strong spin-orbit coupling. The Fermi level is shifted 
upward from the upper slope of j = 5/2 subband in Pu 
and crosses the j = 7/2 subband in Cm corresponding 
to the increasing number of / electrons. A separation of 
the centers of gravity of these subbands for the value of 
Coulomb parameter U — 4 eV results in 5 - 5.5 eV. 



For comparison we present DOSes of the cubic phase 
under pressure calculated within LDA-)-C/+SO for the 
corresponding volumes per actinide atom. As one can 
see, cubic phases with the fitted volumes can be used as 
a good approximation for real phases 5/ DOSes, since 
the centers of gravity and bandwidth of j = 5/2 and 7/2 
subbands are found in good agreement with that in the 
real phases. Also the total density of states of real and 
corresponding cubic phases at the Fermi level are close. 
Having this close similarity in mind, below we report the 
resistivity model results for the DOSes of cubic phases 
with different volumes taken as a starting point. Since 
the model allows to estimate temperature dependence 
of ER, the starting DOSes were temperature broadened. 
The result of such broadening for curium in the volumes 
per atom, corresponding to real volumes are shown in 
Fig. [31 All other DOSes before broadening are reported 
in Ref. 55 and look similar with the broadening. 
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FIG. 3: (Color online) Partial 5/ densities of states of some 
fcc-curium volumes under pressure the CPA broadened for the 
temperatures T = 100, 300, and 500 K. Here the volumes are 
related to Vo - the volume per ion in Cml phase at ambient 
pressure. The Fermi energy corresponds to zero. 



B. Electrical resistivity temperature dependencies 
of metals at normal conditions and under pressure 

Electrical resistivity in CPA is usually calculated 
within Kubo formula for the diagonal part of conduc- 
tivity tensor 



(.r'^i:S&/-(-f)^-^(-) 



(18) 



where 



T{U) = J dEg,{E)E 



7.(^) 



-\ 2 



■(19) 

and gj is DOS of the j-th conductivity band. In this 
equation, the approximate expression v'^ = 2E/m* for 
the square of electron velocity was used. Note, that ap- 
parent limitations of Eq. ()18p arise from neglecting of 
back /(d) — ?> s transitions of conductivity electrons^ 

Numerical solution of Eqs. p^ and ([T5|) with s, d, and 
/ DOSes of metal provides a usual way to estimate effec- 
tive mass of conductivity electrons (m*) from ab initio 
results. Accounting for effective mass renormalization 
and "accepting" DOSes modification with temperature, 
one can try to estimate the absolute values of ER of ac- 
tinides. However, this result will comprise no correction 
on effects of irradiation. 

Full solution of the CPA Eq. (|T5|) for pure fcc-Pu, Am, 
Cm and bcc-Np metals at normal conditions and un- 
der pressure were obtained using corresponding ab initio 
DOSesSS, as the starting point of iteration procedure. 
Bloch constant value equal to 0.8 Ep^ and experimental 



data for the Debye temperatun 



.30.32 



Od -- lOOK, veloc- 



ity of sound and structural data (see Table 1) were used 
for the parametrization of these equations. Note, that 



the Debye temperature estimation a^°i^^'^^ presented pre- 
viously for (5-Pu by other authors differ significantly from 
each other. However, simulations within self-consistent 
general thermodynamic model, accounting for effect of 
anharmonicity of lattice^"^ gives the same value for De- 
bye temperature as found experimentally in Ref. |30. 

It is well known that the value of DOS and its behavior 
in the vicinity of the Fermi level affect significantly ki- 
netic properties of metal. All DOSes modifications with 
temperature and pressure has attracted a great deal of 
attention. The present results are based on the com- 
mon numerical solution of a set of CPA equations (|15l) . 
performed for different pressure and temperature and 
demonstrate general trend of strong influence of electron 
phonon interaction on the initial DOS. In Fig.[3]the den- 
sities of states for a few volumes of curium are shown 
to illustrate the temperature broadening of DOSes. One 
can see from Fig. [31 that strong electron-phonon coupling 
leads to significant smoothing of all initial fine features 
of the DOS curves. At the same time, applied pressure is 
slightly hindered the smoothing of the initial curve but 
does not lead to qualitatively different result. In vicin- 
ity of the melting point, calculated DOSes of all metals 
completely lose all their original features and are similar 
to each other. Calculated dynamics of DOSes evolution 
determines mainly the behavior of ER vs. temperature 
and applied pressure. 

Recently, ER calculations for pure bcc-Np, fcc-Pu, Am, 
and Cm metals within CPA for two band conductivity 
model at high temperatures and normal conditions pro- 
posed typical metallic behavior of ER over the whole tem- 
perature region without any specific peculiarities.^"^ Sim- 
ilar behavior of ER vs. temperature without any anoma- 
lies and singularities were calculated in this work for ac- 
tinides within the proposed conductivity model both at 
normal conditions and under pressure. For bcc-Np and 
fcc-Cm metal only the phonon part of ER was calculated. 
For pure fcc-Pu and Am at high temperature ER is de- 
fined by electron-phonon scattering mainly. 

One can see, that all theoretical curves of ER are sim- 
ilar to the common curves measured for a number of 4(i- 
(5d-) transition metals^ and agree well with the experi- 
mental data. Calculated dependencies of ER for all met- 
als show, that the TCR of the metals increases with the 
pressure and decrease with temperature for 100-^500 K. 
In the vicinity of melting point TCR for all metals was 
found to be a weakly increasing function of temperature. 

Np metal. The first detailed ER experimental data 
on temperature dependence of pure orthorhombic a- 
Np phase with i?295/^4.2=34.36 was reported in 1963.— 
Subsequent experiments (see Ref. [l^ were performed for 
polycrystalline samples of a-Np in 4.2 - 300 K tempera- 
ture range. The residual resistivity of the Np metal was 
found to be equal to 12.2 fiQcm and i?273/^4.2 = 8.15fi^ 
that shows high concentration of different impurities in 
the samples. Authors of Ref. ,20 reported Np R27^/ Ri.2 
value to be equal to 4.47 and pointed out strong effect of 
self-damage on the observed value of ER. 
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FIG. 4: (Color online) Electrical resistivity of Np metal under 
pressure. For each curve the corresponding volume is related 
to the volume per actinide ion (Vo) at ambient pressure - 

V/Vo. 



Note, that residual resistivity values about 
10-^15 fiflcm are typical in concentrated Sd- and 
5(i-transition metal alloys with strong electron impurity 
interaction. After subtraction of this background the 
temperature dependence of ER of Np was found to be 
ordinary for dilute alloys with p '-^ T"^ in low temperature 
region^^ and a weakly increasing non-linear function of 
temperature T > 6d- The calculated ER values agree 
well with experimental data in arbitrary units at normal 
conditions and predicts TCR increasing under pressure. 
Calculations of conductivity of s{p)-, d-, and /-electrons 
show that total ER is determined by s(p) -electrons 
mainly and partially by d-electrons. Assuming that 
/-band electrons are conducting ones, we obtained 
strongly overestimated ER values. Determining the 
temperature-dependent part of ER as p{T) — po and 
using calculated values of ER of s- and d-electrons, we 
found the resistivity values underestimated by 9-11 %. 
That is why only electron-phonon interaction was taken 
into consideration and electron-electron coupling was 
neglected. 

Pu metal. Negative TCR value in a-Pu were mea- 
sured about five decades agoJ^ This result is unique for 
pure metal and was confirmed by several groups J^iii Ex- 
perimental results for ER motivated great discussion of 
possible mechanism of scattering which can provide such 
values of negative coefficient of resistivity and a model 
explaining it. 

However, note at first that a-Pu as well as S-Pu have 
very high values of residual resistivity and high concen- 
trations of impurities and defects. Strong change of 
ER values during the holding time also evidences for 
strong contributions of electron-impurity and electron- 
defect scattering parts to the total value of resistivity. 
Second, thus experimental data were obtained for dilute 
Pu-based alloys and quasi-polycrystalline samples. 

Using well-known experimental fact that ER of poly- 



FIG. 5: (Color online) Electrical resistivity of Pu metal under 
pressure. 



crystalline metals shows behavior analogous to their fcc- 
and bcc-crystal phases, the negative TCR can easily 
be explained in terms of electrons scattering on non- 
coherent phonons. Indeed, electron scattering on lattice 
defects, produced by self-damage in a-Pu and ion mass 
modification as a result of irradiation leads to randomly 
distributed "impurity" ions on the sites of crystal lattice. 
Hence, conductivity electrons scatter not only on "pure" 
phonons but also on randomly distributed Coulomb fields 
of impurities ions. This model equivalence corresponds to 
interference model of scattering in fee- dilute Pu-based al- 
loys ((S-Pu), previously proposed and discussed in Refs.l23l 
and 1 311 . This consideration can be extended, allowing a 
direct way to explain complex experimental data as well 
as non-magnetic nature of the observed ER anomalies, 
p ~ r^ at low temperature, negative TCR at high tem- 
perature range and exclude a-Pu from list of anomalous 
metals, at least concerning its resistivity properties. 

Electrical resistivity calculations of fcc-Pu within pro- 
posed many-band conductivity model were performed in 
absolute units. Experimental ER data for pure fcc-Pu are 
unknown. Numerical simulation shows ordinary temper- 
ature dependence of ER at high temperature. An at- 
tempt to separate resistivity of different groups of elec- 
trons gives domination of s-band electrons with effective 
mass above 25-^30 at high pressure and ~ 35-;-40 at nor- 
mal condition. The electron mass renormalization is the 
result of intensive s — >■ d and s — >■ / transitions and 
strong effect of direct d — > / and back / — > d transitions 
on the DOS values at the Fermi level. Calculated abso- 
lute values of ER underestimate experimental data for 
a-Pu and 6-Pu by a factor of 3 and predict metallic type 
of temperature dependencies of resistivity at normal con- 
ditions and under pressure. TCR of Pu was calculated 
as a weak increasing function of pressure and weak de- 
creasing function of temperature. 

Am metal. During the last five years high interest in 
resistance properties of Am arose due to superconductiv- 
ity found in pure metal. Reported values of Tc of Am 
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FIG. 6: (Color online) Electrical resistivity of Am metal under 
pressure. 



FIG. 7: (Color online) Electrical resistivity of Cm metal under 
pressure. 



metal is above 0.5 K and found to be dependent on ap- 
plied pressure. ^^ 

Electrical resistivity of pur e Am metal was reported 
for the first time in Ref. [la The samples were rela- 
tively pure and low temperature behavior of Am metal 
was obtained as piT) ^ y4±o.5^ that agrees well with 
the data for transition metals containing small number 
of impurities. At high temperature ER demonstrates also 
ordinary dependence with decreasing TCR for all inves- 
tigated temperatures. 

Calculations performed within the proposed model at 
normal conditions and pressure show ordinary ER tem- 
perature dependencies and agree well with experiments, 
see Fig. [S] Obtained TCR is a decreasing function of 
temperature in all regions but slightly increasing with 
applied pressure. Numerical values of ER underestimate 
experimental data by 10 ~ 15 %. Nature of the obtained 
ER dependencies is analogous to calculated for Np and 
Pu ones and determined mainly by renormalization of 
DOSes with temperature. 

At this time nothing can be say about giant values - 
~ 400 iikilcm of ER of Am-IV, obtained at room tem- 
perature and ambient pressure 25-30 GPa.^^ Moreover, 
these giant values of ER were measured in the system 
with typical metal type of conductivity without any sign 
of metal-semiconductor transition. All possible model's 
assumptions, accounting for different scattering mecha- 
nisms of conductivity of electrons and very optimistic 
expectations give ER values about 150-200 fik^cm only. 
The nature of high resistivity state can be explained as- 
suming for example, that heavy electrons of d- and /- 
bands are the conductivity electrons and s-type electrons 
are excluded from charge transport process. However, we 
have no corresponding experimental data to point out 
this opinion. 

Cm metal. Only recently experimental data on ER 
in a-Cm (Cm-I) were reported. Total ER was found as 
a sum of residual resistivity (more than 40 /ifcfi • cm), 
magnetic part of resistivity, and the part of ER origi- 



nated from electron-phonon couplingi^ ER behavior of 
Cm is typical for antiferromagnetic metal and transpar- 
ent anomalies in the vicinity of the Neel point. The mag- 
netic part of ER was determined as a difference between 
the total resistivities of Cm and Am. Note first, that the 
large value of residual resistivity shows that the investi- 
gated samples of Cm metal contained a substantial num- 
ber of impurities, including non-controlled ones. Second, 
nothing can be said about reliability of the suggestion, 
proposed in Ref. [6a, that total ER is a sum of additive 
values, since strong electron-impurity interaction can sig- 
nificantly correct this result. 

In the present paper only the "phonon" part of the 
total ER in Cm was determined and discussed unam- 
biguously. The calculated ER is a usual linear function 
of temperature with positive TCR in the interval of 100 - 
350 K and has similar behavior as in Pu and Am at the 
higher temperature. Temperature coefficient of resistiv- 
ity of Cm demonstrates small decreasing under pressure, 
the similar one was calculated for the other metals un- 
der consideration. For this reason in Cm the DOS value 
at the Fermi level does not change drastically with tem- 
perature and mainly the T/Q factor determines the tem- 
perature behavior of ER. Note, that from 500 to 700 K 
the ER of metals has weak temperature dependence and 
reaches its high temperature limit due to strong erosion 
of DOS at the Fermi level. 



V. CONCLUSION 

We report the results of model calculations of electrical 
resistivity for Np, Pu, Am, and Cm metals under pres- 
sure within CPA for many-bands conductivity model. We 
used DOSes of the metals under consideration as a start- 
ing point in our model calculations. For this purpose, 
Np, Pu, Am, and Cm pure metals were investigated in 
real crystal phases within the LDA-I-C/ method with spin- 
orbit coupling (LDA+f7+S0). For Am and Pu metals in 



all structures we found nonmagnetic ground state. We 
compared DOSes calculated for the real crystal struc- 
tures of the transuranium metals with the DOSes for 
corresponding fee volumes per actinide atom. It follows 
from this comparison that total density of states, as well 
as centers of gravity and bandwidth of j — 5/2 and 7/2 
subbands are in good agreement. It means that the cubic 
structures in different volumes provide good estimation 
for the DOSes of real crystal structures. 

Results of our CPA calculations within derived many- 
band conductivity model for pure fcc-Pu, Am, and Cm 
and bcc-Np using Bloch constant value ^ O.SEp^'^ and 
ab initio DOSes of metals show ordinary metal type of 
ER behavior vs. temperature at normal conditions and 
under pressure. High values of ER in these metals are a 
consequence of s — > d and s — )■ / transitions of scattering 
electrons and effective mass of conductivity electrons in- 
creasing as a result of s- and d-electron bands hybridiza- 
tion. Weak non-linearity in the ER temperature behavior 
of metals above 9d is caused by DOS erosion at the Fermi 
level due to electron-phonon interaction. 

The obtained high-resistivity values in actinides result 
from s — >■ d and s — > / interband transitions and conduc- 
tivity electron mass renormalization due to strong s-d hy- 
bridization. Also strong influence on the ER temperature 
dependencies of actinides is a result of DOS of accepting 
bands modifications as products of direct s — >■ d and 
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d — >■ / and back / — >■ d electrons transitions. Calculated 
phonon part of resistivity underestimates experimental 
data by a factor of 1.5 for Cm and by a factor of 3 for 
ideal Pu. 
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VII. APPENDIX 

Let us assume that the matrix elements of interaction 
Vw are averaged over the angle between the wave vectors 
k and k' and thus depend only on the band index. Using 
this simplification and Eq. (|13p. one obtains: 



1 






J(k-k',R„) 



Y.i^jjAumAV,,riu)] ++J2[V,,Au)]F,,[V,,,,iu)]F,,[V,,,,iu)] + ■■■} Swa+a,, (20) 



In I Ji Jlj2 

where the j'-th band electron Green function Fj is defined as 



^J - w z^ ■ 



k '•' 

Comparing Eqs. (|S]) and (|20p . for Aj one obtains 



(21) 



(22) 



Using the following matrix form 



[F] 



"s 1 




Fd 


, [A] = 


OOF/ 





A,, 





Arf 








A/ J 



, [V] 



Vn^ss{u) Vn,sdiu) Vn,sf{u) 
Vn,dsiu) Vn,ddiu) Vn,df{u) 
Vnjs{u) Vnjdiu) Vnjf{u) 



(23) 



the series (IT^ is summed up accurately in the convergence range |Vj(m)Fj| < 1. Within single-electron and single-site 
approaches that gives 



m) = {[V][F][V] + [V][F][V][F][V] 



1-[F][V]/- 



(24) 



Calculating matrix products and averaging over phonons both part of Eq. (|24|) . the set of CPA equations can be 
obtained. 



11 



* Electronic address: 'y.tsiovkin@mail.ustu.ru' 

^ Challenges in Plutonium Science, edited by N. G. Cooper, 
Los Alamos Sci. 26 (LANL, Los Alamos, NM, 2000). 

^ D. A. Young, Phase diagram of the Elements (University 
of California Press, Berkley, 1991). 

^ S. S. Hecker, D. R. Harbur, and T. G. Zocco, Prog. Mater. 
Sci. 49, 429 (2004). 

* G. H. Lander, Science 301, 1057 (2003). 

^ S. Heathman, R. G. Haire, T. Le Bihan, A. Lindbaum, 
M. Idiri, P. Normile, S. Li, R. Ahuja, B. Johansson, and 
G. H. Lander, Science 309, 110 (2005). 
® S. Heathman, R. G. Haire, T. Le Bihan, A. Lindbaum, 
K. Litfin, Y. Meresse, and H. Libotte, Phys. Rev. Lett. 
85, 2961 (2000). 
'^ A. Lindbaum, S. Heathman, K. Litfin, Y. Meresse, 
R. G. Haire, T. Le Bihan, and H. Libotte, Phys. Rev. B 
63, 214101 (2001). 
® G. van der Laan, K. T. Moore, J. G. Tobin, B. W. Chung, 
M. A. Wall, and A. J. Schwartz, Phys. Rev. Lett. 93, 
097401 (2004); G. van der Laan and B. T. Thole, Phys. 
Rev. B 53, 14 458 (1996). 
^ K. T. Moore, M. A. Wall, A. J. Schwartz, B. W. Chung, 
D. K. Shuh, R. K. Schulze, and J. G. Tobin, Phys. 
Rev. Lett. 90, 196404 (2003); K. T. Moore, M. A. Wall, 
A. J. Schwartz, B. W. Chung, S. A. Morton, J. G. Tobin, 
S. Lazar, F. D. Tichelaar, H. W. Zandbergen, P. Soderlind, 
and G. van der Laan, Philos. Mag. 84, 1039 (2004). 

^° K. T. Moore and G. van der Laan, Rev. Mod. Phys. 81, 
235 (2009). 

" K. T. Moore, G. van der Laan, R. G. Haire, M. A. Wall, 
A. J. Schwartz, and P. Soderlind, Phys. Rev. Lett. 98, 
236402 (2007). 

^^ H. Ledbetter, A. Lawson, and A. Migliori, J. Phys.: Con- 
dens. Matter 22, 165401 (2010). 

" A. M. Boring and J. L. Smith, Los Alamos Sci. 26, 90 
(2000); S. S. Hecker and J. C. Martz, ibid. 26, 238 (2000). 

" R. Smoluchowskii, Phys. Rev. 125, A1577 (1962). 

^^ C. E. Olsen and R. O. Elliott, Phys. Rev. 139, A437 
(1965). 

1® M. B. Brodsky, Phys. Rev. 137, A1423 (1965). 

^■^ M. B. Brodsky, Rep. Prog. Phys. 41, 1547 (1978). 

^* W. MiiUer, R. Schenkel, H. E. Schmidt, J. C. Spirlet, 

D. L. McElroy, R. O. A. Hall, and M. J. Mortimer, J. Low 
Temp. Phys. 30, 561 (1978). 

^^ G. T. Meaden, Proc. R. Soc. Lond. A 276, 553 (1963). 

^° E. King, J. A. Lee, K. Mendelssohn, and D. A. Wigley, 
Proc. R. Soc. Lond. A 284, 325 (1965). 

^^ It is a well-known fact, that the order in p ~ T" dependen- 
cies increases up to n = 5 with purity of transition metals. 

^^ J.-C. Griveau, J. Rebizant, G. H. Lander, and G. Kotliar, 
Phys. Rev. Lett. 94, 097002 (2005). 

^^ Yu. Yu. Tsiovkin, M. A. Korotin, A. O. Shorikov, 
V. I. Anisimov, A. N. Voloshinskii, A. V. Lukoyanov, 

E. S. Koneva, A. A. Povzner, and M. A. Surin, Phys. Rev. 
B 76, 075119 (2007). 

^^ S. Meot-Reymond and J. M. Fournier, J. Alloys Compd. 

232, 119 (1996). 
^^ V. Dallacasa, J. Phys.: F Metal Phys. 11, 177 (1981). 
^^ M. Fluss, B. D. Wirth, M. Wall, T. E. Felter, M. J. Caturla, 

A. Kubota, and T. Diaz de la Rubia, J. Alloys Compd. 

368, 62 (2004). 



^■^ R. JuUien, M. T. Beal-Monod, and B. Cogblin. Phys. Rev. 

B 4, 1441 (1974). 
'^^ Yu. Piskunov, K. Mikhalev, A. Gerashenko, A. Pogudin, 

V. Ogloblichev, S. Verkhovskii, A. Tankeyev, V. Arkhipov, 

Yu. Zouev, and S. Lekomtsev, Phys. Rev. B 71, 174410 

(2005). 
^^ S. V. Verkhovsku, V. E. Arkhipov, Yu. N. Zuev, 

Yu. V. Piskunov, K. N. Mikhalev, A. V. Korolev, I. L. Svy- 

atov, A. V. Pogudin, V. V. Ogloblichev, and A. L. Bu- 

zlukov, JETP Lett. 82, 139 (2005). 
^° J. C. Lashley, J. Singleton, A. Migliori, J. B. Betts, 

R. A. Fisher, J. L. Smith, and R. J. McQueeney, Phys. 

Rev. Lett. 91, 205901 (2003). 
^^ Yu. Yu. Tsiovkin and L. Yu. Tsiovkina, J. Phys.: Condens. 

Matter 19, 056207 (2007). 
^^ J. C. Lashley, A. C. Lawson, R. J. McQueeney, and 

G. H. Lander, Phys. Rev. B 72, 054416 (2005). 
^^ A. O. Shorikov, A. V. Lukoyanov, M. A. Korotin, and 

V. I. Anisimov, Phys. Rev. B 72, 024458 (2005). 
^^ A. B. Shick, V. Drchal, and L. Havela, Europhys. Lett. 69, 

588 (2005). 
'^^ A. B. Shick, J. Kolorenc, A. L Lichtenstein, and L. Havela, 

Phys. Rev. B 80, 085106 (2009). 
^^ L. Havela, A. Shick, and T. Gouder, J. Appl. Phys. 105, 

07E130 (2009). 
^■^ R. Atta-Fynn and A. K. Ray, Europhys. Lett. 85, 27008 

(2009). 
^^ F. Cricchio, F. Bultmark, and L. Nordstrom, Phys. Rev. 

B 78, 100404 (2008). 
'^^ F. Bultmark, F. Cricchio, O. Granas, and L. Nordstrom, 

Phys. Rev. B 80, 035121 (2009). 
*° G. Kotliar, S. Y. Savrasov, K. Haule, V. S. Oudovenko, 

O. ParcoUet, and C. A. Marianetti, Rev. Mod. Phys. 78, 

865 (2006). 
*^ S. Y. Savrasov, G. Kotliar, and E. Abrahams, Nature (Lon- 
don) 410, 793 (2001). 
*^ L. V. Pourovskii, M. I. Katsnelson, A. L Lichtenstein, 

L. Havela, T. Gouder, F. Wastin, A. B. Shick, V. Drchal, 

and G. H. Lander, Europhys. Lett. 74, 479 (2006). 
*^ L. V. Pourovskii, G. Kotliar, M. L Katsnelson, and 

A. I. Lichtenstein, Phys. Rev. B 75, 235107 (2007). 
*'' J. H. Shim, K. Haule, and G. Kotliar, Nature (London) 

446, 513 (2007). 
^^ J. H. Shim, K. Haule, and G. Kotliar, Europhys. Lett. 85, 

17007 (2009). 
""^ J.-X. Zhu, A. K. McMahan, M. D. Jones, T. Durakiewicz, 

J. J. Joyce, J. M. Wills, and R. C. Albers, Phys. Rev. B 

76, 245118 (2007). 
"■^ C. A. Marianetti, K. Haule, G. Kothar, and M. J. Fluss, 

Phys. Rev. Lett. 101, 056403 (2008). 
"^ J. G. Tobin, P. Soderlind, A. Landa, K. T. Moore, 

A. J. Schwartz, B. W. Chung, M. A. Wall, J. M. Wills, 

R. G. Haire, and A. L. Kutepov, J. Phys.: Condens. Mat- 
ter 20, 125204 (2008). 
**' V. I. Anisimov, F. Aryasetiawan, and A. I. Lichtenstein, 

J. Phys.: Condens. Matter 9, 767 (1997). 
^° V. I. Anisimov and O. Gunnarsson, Phys. Rev. B 43, 7570 

(1991). 
^^ O. Gunnarsson, O. K. Andersen, O. Jepsen, and J. Zaanen, 

Phys. Rev. B 39, 1708 (1989). 
^^ All calculations were performed in the TB-LMTO-ASA 



12 



package. O. K. Andersen, Phys. Rev. B 12, 3060 (1975). 

^^ A. N. Voloshinskii and A. G. Obukhov, Phys. Met. Met- 
tallogr. 91, 238 (2001). 

^"^ A. B. Chen, G. Weisz, and A. Sher, Phys. Rev. B 5, 2897 
(1972). 

^^ A. V. Lukoyanov, A. O. Shorikov, V. B. Bystrushkin, 
A. A. Dyachenko, L. R. Kabirova, Yu. Yu. Tsiovkin, 
A. A. Povzner, V. V. Dreniov, M. A. Korotin, and 
V. I. Anisimov, larXiV: 1005.2851. 

^® For the phases with several inequivalent actinide ions we 
found close values of the off-diagonal elements and mo- 
ments. For such phases the averaged values per ion and 
the largest OD elements are given. 

^"^ W. H. Zachariasen, Acta Crystallogr. 5, 660 (1952). 

^^ W. H. Zachariasen, Acta Crystallogr. 5, 664 (1952). 

^^ W. B. Pearson, Handbook of Lattice Spacings and Struc- 
tures of Metals and Alloys, (Pergamon, New York, 1967), 



Vol. 2. 
®° W. H. Zachariasen and F. H. EUinger, Acta Crystallogr. 

8, 431 (1955); ibid. 16, 777 (1963). 
^^ P. G. Huray and S. E. Nave, in Handbooks on the Physics 

and Chemistry of the Actinides, edited by A. J. Freeman 

and G. H. Lander (Elsevier, Amsterdam, 1985), Vol. 5, 

p. 311. 
"^ H. Ledbetter, A. Migliori, J. Belts, S. Harrington, and 

S. El-Khatib, Phys. Rev. B 71, 172101 (2005). 
^'^ A. N. Filanovich, A. A. Povzner, V. Yu. Bodryakov, 

Yu. Yu. Tsiovkin, and V. V. Dremov, Technical Physics 

Letters 35, 929 (2009). 
^^ T. Aisaka and M. J. Shimizu, J. Phys. Soc. Jap. 28, 646 

(1970). 
®^ The dependence T^-^^ is more correct.— 
^^ R. Schencel, Solid State Comm. 23, 389 (1977). 



